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PREFACE 


This  report  describes  a computer  model  developed  at  the 
Transportation  Systems  Center,  Cambridge,  Massachusetts,  for 
the  simulation  of  electrical  power  distribution  characteristics 
of  Automated  Guideway  Transit  (AGT)  systems.  This  work  was 
conducted  under  the  sponsorship  of  the  Office  of  New  Systems 
and  Automation  of  the  Urban  Mass  Transportation  Administration 
(UMTA) . The  objective  of  this  simulation  effort  is  to  provide 
the  necessary  software  for  rapid  evaluation  and  assessment  of 
AGT  power  system  distribution  requirements. 
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SUMMARY 


This  report  describes  a simulation  model  developed  to 
study  the  electrical  power  distribution  characteristics  of 
Automated  Guideway  Transit  (AGT)  systems.  This  work  was 
conducted  under  sponsorship  of  the  Advanced  Group  Rapid 
Transit  (AGRT)  Program  of  the  Office  of  New  Systems  and 
Automation,  Urban  Mass  Transportation  Administration  (UMTA) . 

The  purpose  of  the  AGRT  program  is  the  development,  evaluation, 
and  verification  of  a second  generation  AGT  technology  which 
is  capable  of  providing  transit  service  in  urban  areas.  The 
computer  simulation  model  developed  at  the  Transportation 
Systems  Center  provides  the  required  software  for  rapid 
evaluation  of  the  complex  power  distribution  of  multi -vehicle 
AGT  systems.  The  application  of  the  computer  model  provides 
the  system  designers  with  important  information  on  the 
anticipated  power  consumption  characteristics  of  AGT  systems. 
Such  information  is  extremely  valuable  in  both  the  initial 
design  stage  of  system  development  as  well  as  in  later  system 
evaluation  studies  aimed  at  assessing  the  overall  system 
performance . 

The  purpose  of  this  report  is  to  provide  the  reader  with 
a working  knowledge  of  the  multi -vehicle  AGT  simulation  model 
and  its  application  to  urban  transit  systems.  Basically  the 
model  is  comprised  of  two  separate  elements  which  describe 
(1)  the  time-varying  spatial  characteristics  of  the  power 
distribution  network,  and  (2)  the  complex  power 

viii 


characteristics  of  the  active  (vehicle)  loads.  These  elements 
are  integrated  into  a comprehensive  simulation  model  which 
accepts  as  prescribed  input  data  the  vehicle  mission  profile 
and  computes  as  output  data  (in  graphical  or  tabular  format) 
the  power  consumption  characteristics  of  the  vehicle  system 
network.  The  report  includes  a complete  listing  of  the 
FORTRAN  program  and  provides  an  illustrative  example  of  its 
application  to  multi-vehicle  ACT  systems. 
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1. 


INTRODUCTION 


This  document  describes  a FORTRAN  computer  program 
which  models  the  electric  power  requirements  of  a guideway 
transit  system.  Typical  systems  which  could  be  modeled 
include:  Morgantown  Personal  Rapid  Transit  System,  Dallas- 

Fort  Worth  Airtrans , or  one  of  the  proposed  Downtown  People- 
Movers.  All  of  these  are  Automated  Guideway  Transit  (AGT) 
systems.  The  model  is  designed  for  studies  which  examine 
the  impact  on  power  consumption  of  the  following: 

1)  Changing  the  vehicle  propulsion  system 

2)  Changing  the  power  distribution  system  by  moving/ 
adding/deleting  power  substations 

3)  Incorporating  regenerative  braking 

4)  Altering  the  vehicles'  mission  profiles. 

The  model  is  general  enough  so  that  different  AGT  systems  can 
be  studied.  The  model  is  a modular  one  so  that  AGT  system 
modifications  are  made  either  by  changing  modules  (FORTRAN 
subroutines)  or  by  changing  parameters  within  a module  (FORTRAN 
statements) . 

The  important  features  of  this  model  are:  (1)  the  equiva- 

lent circuit  for  the  electric  power  distribution  system  is 
automatically  updated  as  the  vehicles  move,  (2)  the  power  flow 
problem  is  solved  using  the  efficient  Newton-Raphson  algorithm, 
(3)  harmonic  currents  (which  are  important  in  sizing  power 
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distribution  equipment)  are  included. 

The  remainder  o£  this  document  is  structured  to  serve  as 
both  a users’  manual  and  a programmers'  reference  manual. 
Section  2 is  an  overview  of  the  model  and  it  contains  a func- 
tional description  of  the  model  and  its  component  parts. 
Section  3 contains  the  analysis  used  in  the  model.  The  simu- 
lation output  is  discussed  in  Section  4.  Section  5 shows  the 
conclusions  of  a model  validation  study  using  the  AIRTRANS 
system . 

Finally  the  procedure  for  preparing  input  data,  and  using 
the  model  are  included  in  Section  6. 
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2.  MODEL  OVERVIEW 


The  power  flow  in  an  ACT  system  is  a complex  process. 

In  a simplified  form,  three-phase  electric  power  enters  the 
ACT  system  via  an  electric  utility  connection.  It  flows 
through  cables  to  power  substations  located  along  the  guide- 
way. At  the  substation,  a step-down  transformer  converts  it 
to  a lower  voltage,  current  is  fed  to  the  power  rail  on 
the  guideway,  and  current  may  be  fed  to  other  buildings 
for  housekeeping  use  (heating,  cooling,  lighting,  etc.) 

The  power  rail  carries  the  current  to  the  AGT  vehicles 
on  the  guideway.  Once  the  electric  power  reaches  the  vehicles, 
most  of  it  is  converted  into  thrust  by  the  propulsion  system, 
but  some  is  used  for  vehicle  housekeeping. 

The  thrust  produced  by  the  propulsion  system  is  a 
function  of  time.  The  vehicle’s  control  system  will  vary 
the  thrust  as  it  adjusts  the  velocity  of  the  vehicle.  As  a 
result  each  vehicle's  power  demand  varies  with  time,  and  the 
AGT  system’s  power  demand  may  also  vary. 

In  effect  the  AGT  is  a power  distribution  network  with 
(1)  generators  (the  electric  utility  connection)  , (2)  trans- 
mission lines  (cables  and  power  rail) , and  (3)  loads  (the  AGT 
vehicles) . 

The  solution  to  the  power  demand  problem  is  to  model  the 
AGT  system  as  a power  distribution  network  and  solve  the  load- 
flow  problem  repeatedly  as  the  time-varying  loads  change  and 
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as  the  transmission  lines  change.  One  point  that  needs  to  be 
made  is  that  in  simplifying  the  ACT  system,  we  neglected 
regenerative  braking.  Regenerative  braking  means  that  during 
braking  operations,  the  ACT  vehicles  can  generate  power  and 
feed  it  back  into  the  power  distribution  system.  Thus  our 
power  distribution  network  has  loads  which  vary  in  magnitude 
and  sign  (negative  loads  are  generators) . A block  diagram  of 
the  ACT  power  distribution  model  is  shown  in  Figure  2-1. 

Note  that  the  three  quantities:  (1)  the  vehicle's  position, 

(2)  the  vehicle's  power  demand,  and  (3)  the  housekeeping  power 
demand,  are  the  independent  variables.  Together  they  determine 
the  systems  aggregate  power  demand. 

Each  vehicle  in  the  AGT  fleet  is  following  a route  along 
the  guideway.  The  characteristics  of  this  route  (i.e.  the 
grades,  the  expected  vehicle  velocity,  and  average  headwinds 
encountered)  and  the  characteristics  of  the  vehicle's  propul- 
sion system  (thrust  or  tractive  force  at  the  wheels  vs.  electric 
power  into  the  motor)  both  determine  the  vehicle's  time-varying 
power  demand.  Figure  2-2  shows  a block  diagram  of  the  AGT 
vehicle  model. 

The  inputs  to  the  vehicle  model,  which  relate  to  the  guideway 
route,  are  functions  of  time.  This  set  of  functions  (expected 
velocity,  guideway  grade,  and  headwind  encountered)  form  part 
of  a vehicle's  mission  profile. 
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FIGURE  2-1.  FUNCTIONAL  VIEW  OF  THE  POWER  DISTRIBUTION  MODEL 
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So  far  the  propulsion  system  model  has  described  how 


much  power  each  vehicle  is  demanding,  but  not  where  in  the 
power  distribution  network  the  demand  is.  After  all,  some 
of  the  ACT  vehicles  are  moving  along  the  power  rail  some- 
where between  stations.  So  we  augment  the  mission  profile 
to  include  a location  code  as  a function  of  time.  The  code 
has  two  parts:  (1)  a prefix  part  to  identify  a guideway 

segment,  and  (2)  a footage  part  to  identify  the  distance  from 
the  segment  reference  end  point. 

By  using  the  location  code  for  vehicles  and  power  sub- 
stations, together  with  information  on  the  spatial  arrange- 
ment of  guideway  segments  into  the  guideway  network,  the 
physical  length  of  the  power  rail  which  separates  vehicles 
and  power  substations  can  be  found.  (See  Figure  2-3  for  an 
example.)  This  physical  length  can  be  converted  into  the  power 
rail's  transmission  line  admittance  for  use  in  the  power  dis- 
tribution model. 

The  transformation  from  substation  and  vehicle  location 
codes  into  power-rail  admittances  is  summarized  by  the  guide- 
way network  model.  That  model  is  shown  functionally  in  Figure 

4. 

In  summary,  the  model  has  three  parts: 

1.  ACT  vehicle  propulsion  system  model--it  calculates 
the  individual  vehicle's  power  demand  versus  time. 
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2,  Guideway  Network  Model  - It  calculates  the  position 
of  the  individual  vehicle  versus  time. 

3.  Power  Distribution  Model  - It  combines  the  individual 

power  demands,  based  on  their  magnitude  and  their 

location,  into  a system  power  demand  versus  time. 

2.1  OUTPUT 

The  model  calculates  a number  of  values  in  addition  to 
the  systems'  electric  power  demand  at  each  sample  time.  In 
general,  the  output  data  will  depend  on  whether  the  power 
rail  carries  ac  or  dc . A general  list  of  data  follows: 

1.  Active  power  (P)  at  the  power  rail/vehicle  inter- 

face. P,  which  is  either  dc  or  60  Hz  ac , is  the  real 
power  used  for  the  vehicles'  propulsion  system  and 
its  auxiliary  loads  (i.e.  heating,  lighting,,  cooling 
and  control  systems) . 

2.  Reactive  power  (Q)  at  the  power  rail/vehicle  inter- 
face. Q,  which  is  present  in  the  ac  power  rail  case, 
is  the  quadrature  component  of  the  60  Hz  real  power. 

3.  The  power  rail  voltage.  This  is  the  voltage  present 
at  the  power  rail/vehicle  interface.  It  is  either 
dc  or  60  Hz  ac . 

4.  Distortion  current  (10)  at  the  vehicle.  It  is  the 
RMS  value  of  the  harmonic  currents  which  flow  across 
the  power  rail/vehicle  interface. 
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5.  Distortion  current  at  the  utility.  The  RMS 
value  o£  the  harmonic  currents  which  flow  from 
the  utility  to  the  ACT  system. 

6.  The  active  power  (P)  and  reactive  power  (Q)  at 
the  utility. 

Assuming  that  the  utility  is  a "stiff”  source  and  therefore  its 
voltage  (El)  is  sinusoidal,  then  the  apparent  power  (U)  at 
the  utility  is  given  by  Equation  2.1-1  (where  **  indicates 
exponentiation) . 

U = SQRT  (p  **  2 + Q **  2 + El  **  2 * ID  **  2)  (2.1-1) 

The  above  values  are  summarized  in  a table  printed  by  the 
model  at  each  sample  time. 

2.2  INPUT 

The  simulation  involves  data  from  two  input  files.  The 
first  is  the  guideway  deployment  data  file.  It  contains  infor- 
mation about  the  guideway  which  is  relevant  to  the  power 
distribution  circuit.  Each  record  in  the  file  describes  a 
branch  in  the  power  distribution  circuit.  The  following  list 
shows  the  type  of  information  in  the  file: 

1.  Power  feed  point  locations 

2.  Power  cross-under  locations.  These  are  the 
points  at  which  the  two  power  rails  are  connected 
together  at  a merge  or  diverge  structures. 

3.  Length  of  power  rail  segments  between  feed  point 
locations  and  cross-under  locations.  These  are  used 
to  calculate  the  impedance  of  the  branches  in  the 
power  distribution  circuit. 
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4. 


Step-down  transformer  leakage  impedance.  If  the 
branch  between  the  utility  and  power  feed  point 
contains  a step-down  transformer,  then  its  impedance 
plus  any  cable  impedance  may  be  included. 

In  summary,  the  guideway  deployment  data  file  describes 
the  nodes  and  branches  (in  the  power  distribution  circuit) 
which  are  geographically  fixed. 

The  second  is  the  mission  profile  file.  It  contains 
information  about  the  vehicles ' locations  and  operating 
conditions.  Each  record  describes  a point  on  an  ideal  mission 
profile.  The  following  list  shows  the  type  of  information  in 
the  file: 

1.  A time  increment.  This  is  the  time  interval 
which  normally  ellapses  before  the  next  mission 
profile  point  describes  the  vehicle  (i.e.,  station 
dwell  time) . 

2.  A commanded  velocity.  This  is  the  velocity 

which  is  an  ideal  velocity,  and  the  vehicle  acceler- 
ates or  decelerates  towards  it. 

3.  Encountered  grade.  This  measures  the  rise  or 
descent  in  the  guideway  encountered  by  the  vehicle. 
Note  that  if  superelevation  is  present,  then  it  can 
be  combined  with  the  grade  to  produce  a compensated 
grade . 

4.  Encountered  headwind. 
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5,  Vehicle  location.  This  is  the  position  of  the 

vehicle  on  the  guideway  where  the  mission  profile 
point  data  becomes  effective. 

It  should  be  noted  that  the  mission  profile  is  an  ideal 
one  from  which  a vehicle  may  deviate.  For  example,  if  the 
vehicle's  propulsion  system  cannot  provide  the  necessary 
tractive  effort  or  if  it  is  j erk- limited , then  the  current 
velocity  differs  from  the  commanded  velocity. 

2.3  SUBMODELS 

This  section  considers  the  three  submodels  and  their 
interfaces.  The  model  can  be  considered  a framework  into 
which  modules  are  inserted.  The  user  creates  these  modules 
for  a particular  vehicle  and  a particular  guideway  deployment. 
Figure  2-5  shows  the  three  models  and  their  interfaces. 

A functional  description  for  each  submodel  follows. 

2.3.1  Vehicle  Propulsion  Submodel 

This  model  calculates  the  electrical  power  load  which 
the  vehicle  produces  at  the  power  rail.  This  load  is  a function 
of  the  vehicle  dynamics,  motor,  the  power  conditioning  unit, 
and  the  auxiliary  loads.  The  computation  of  the  load  rises 
from  the  following  five  steps: 

1.  Data  on  a vehicle's  commanded  mission  profile 

(velocity,  acceleration,  encountered  grade)  and  data 
on  the  vehicle's  current  state  (actual  velocity  and 
position)  are  retrieved  from  tables. 
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FIGURE  2-5.  THREE  SUBMODELS  AND  THEIR  INTERFACES 


2.  The  equations  o£  motion  for  the  vehicle  (which 
include  a quadratic  polynominal  in  velocity  to 
model  tractive  resistance)  are  solved  for  the 
mechanical  power  of  the  motor. 

3.  The  equations  for  the  motor  are  solved  for  the 
electric  power  into  the  motor  terminals. 

4.  The  equations  for  the  power  conditioning  unit  are 
solved  for  the  electric  power  load  at  the  power 
rail.  These  equations  use  as  input  the  power  at 
the  motor  terminals. 

5.  Data  on  the  vehicles'  load  at  the  power  rail/ 
vehicle  interface  and  at  the  motor  terminals  is 
entered  into  tables  which  are  accessible  to  the 
other  submodels. 

It  should  be  noted  that  steps  1 and  5 are  bookkeeping 
operations  required  to  keep  track  of  a fleet  of  vehicles.  The 
user  is  only  concerned  with  providing  information  about  the 
dynamics,  the  motor  system,  and  the  power  conditioning  system 
of  a single  vehicle.  The  computations  are  automatically 
extended  to  the  fleet  of  vehicles  and  tables  are  used  to  store 
the  fleet  data. 

2.3.2  Guideway  Deployment  Submodel 

This  model  is  used  for  bookeeping.  In  general,  it  is  a 
table  of  power  rail  segments.  The  data  recorded  for  each 
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segment  includes  its  length  and  its  adjoining  segments.  The 
computation  done  by  this  submodel  entails  updating  a queue  o£ 
the  vehicles  on  each  power  rail  segment.  This  is  done  using 
the  current  position  of  each  vehicle. 

2.3.3  Power  Flow  Submodel 

This  model  solves  the  equivalent  circuit  for  the  power 
distribution  network.  It  is  an  iterative  computation,  similar 
to  the  algorithm  used  by  electric  utility  companies  to  simulate 
their  power  transmission  networks.  The  difference  here  is  that 
the  loads  are  dynamic.  That  is,  they  vary  with  the  voltage 
at  the  vehicle/power  rail  interface.  The  power  flow  submodel 
involves  the  following  computational  steps:  (in  this  model,  we 

use  the  name  "node”  where  nodes  are  power  feed  points,  power 
rail  connections,  or  vehicle/power  rail  interface  points). 

1.  Using  the  data  from  the  guideway  deployment  table 
(which  also  specifies  the  vehicles  located  on  each 
segment  at  a given  time) , calculate  the  node 
admittance  matrix. 

2.  Assume  initial  node  voltage  values.  Calculate  the 
real  and  reactive  power  at  the  vehicle/power  rail 
interface.  This  is  the  power  needed  at  each  vehicle. 
Call  it  the  scheduled  power. 

3.  Set  an  iteration  counter  to  zero. 

4.  Calculate  the  power  flowing  in  the  power  rail.  From 
this  the  power  delivered  at  each  vehicle  is  determined. 
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5. 


Find  the  largest  discrepancy  between  the  power 
scheduled  and  power  delivered. 

6.  I£  the  discrepancy  is  small  enough  (less  than  one 
percent,  say),  then  end  the  process  with  step  7. 
Otherwise  go  on  to  step  8. 

7.  Print-out  the  voltages  and  power  values.  The  power 
distribution  computation  ends  here. 

8.  Calculate  a Jacobian  matrix  (its  elements  are  partial 
derivatives  o£  power  with  respect  to  voltage) . 

9.  Invert  the  Jacobian  matrix. 

10.  Use  it  to  calculate  the  updated  node  voltages. 

11.  Recalculate  the  scheduled  power  at  the  power  rail/ 
vehicle  inter£aces.  Continue  at  step  4. 

In  Section  3,  the  details  and  analysis  used  £or  these  models 
are  presented. 


3.  ANALYSIS  FOR  THE  MODELS 


In  this  section,  the  data  and  analysis  for  a test  case 
are  discussed.  The  program  listings  in  Appendix  A have  been 
set  up  for  this  same  test  case.  In  general,  the  vehicle  is 
externally  similar  to  the  Morgantown  vehicle,  but  with  a 120 
horsepower  dc  motor  supplied  by  three-phase  power  rail 
through  a variable-voltage  rectifier.  The  guideway  is  a loop 
which  is  based  on  the  Morgantown  system.  Two  vehicles  are 
moving  around  the  loop  following  a mission  profile  which 
approximates  the  Morgantown  operation  (but  with  only  two 
vehicles  out  of  the  usual  fifteen  vehicle  fleet  running) . The 
analysis  is  divided  into  seven  parts  - one  for  each  major 
subroutine . 

3.1  VEHICLE  DYNAMICS 

In  this  section  the  dissipative  and  conserved  forces 
acting  on  the  vehicle  are  used  to  calculate  the  real  power 
at  the  wheels.  When  positive,  this  power  must  be  supplied  by 
the  propulsion  system.  When  negative,  this  power  must  be 
absorbed  by  friction  brakes  or  some  power  source  (regenerative 
braking) . 

Rolling  friction  force  arises  because  of  the  roughness 
of  the  contact  between  the  moving  vehicle  and  the  guideway. 

It  is  linearly  proportional  to  the  vehicle  weight  with  an 
experimentally  determined  proportionality  constant. 
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VM  = 12,000  vehicle  weight  (lb) 

A = 0.025  coefficient  of  rolling  friction  (n.d.) 

GRADE  = encountered  grade  in  percent 

FR  = VM  * COS  (GRADE/100.)  * A where  FR  is  the  rolling 

friction  force  (lb) 

Coulombic  friction  forces  arise  because  of  wheel/guideway 
deformation  and  wheel  bearing  resistance.  It  is  proportional 
to  the  vehicle  weight  and  speed. 

VV  = vehicle  velocity  (mph) 

VM  = 12,000  vehicle  weight  (lb) 

B = 0.00005  coulombic  friction  coefficient  (n.d.) 

FC  = VM  * VV  * B where  FC  is  the  coulombic  friction  force  (lb) 
Aerodynamic  drag  forces  arise  because  of  vehicle  motion 
relative  to  the  air  mass.  It  is  proportional  to  frontal  area 
and  the  vehicle's  relative  velocity  squared. 

C = 0.85  drag  coefficient  for  leading  vehicle,  use  0.19 
for  trailing  vehicles 

RHO  = 0.002331  air  density  (slugs/ft  **  3) 

CONI  = 1.46667  ft/sec  per  mph,  conversion  factor 
W = vehicle  velocity  (mph) 

HW  = encountered  headwind  (mph) 

FD  = 0.5  * RHO  * C * (VV+HW)  **  2*  CON  1 **  2 
where  FD  is  aerodynamic  drag  force  (lb) 

The  total  dissipative  force  can  be  found  as  RFORCE. 

RFORCE  = FR  + FC  + FD 

The  gravitational  force  arises  because  of  a component  of 
the  vehicle's  weight  retards  its  uphill  motion  and  advances  its 
downhill  motion.  It  is  calculated  as  the  component  of  weight 
tangent  to  the  guideway. 
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VM  = 12,000  vehicle  weight  (lb) 

GRADE  = encountered  grade 

EG  = VM  * SIN  (ATAN  (GRADE/100))  where  EG  is  the 
gravitational  force  (lb) 

Inertial  force  arises  when  the  vehicle  must  change  velocity 
and  overcome  inertia.  It  is  proportional  to  the  vehicle's 
mass  and  acceleration. 

VM  = 12,000  vehicle  weight  (lb)l 
KG  = 32.174  gravity  (Ib/slug) 

ACC  = vehicle  acceleration  (mph/sec) 

CONI  = 1.46667  ft/sec  per  mph  Conversion 

El  = VM/G  * ACC  * CONI  where  El  is  inertial  force  (lbs) 

The  total  tractive  force,  ET,  can  be  calculated.  It  is 
the  thrust  required  at  the  wheels  in  order  to  maintain  a given 
velocity  and  acceleration  with  an  encountered  guideway  grade 
and  headwind. 

ET  = REORCE  + EG  + El  where  ET  is  the  tractive  force  (lb) 
The  power,  PW,  required  to  develop  the  necessary  tractive 
force  at  the  vehicle  wheels  is  the  product  of  that  force  and 
vehicle  velocity. 

ET  = tractive  force  at  wheels  (lb) 

W = vehicle  velocity  (mph) 

CONI  = 1.46667  ft/sec  per  mph  Conversion  factor 
C0N2  = 1.356  watt  per  ft  Ib/sec  Conversion  factor 
PW  = PT*W*C0N1  CON2  where  PW  is  the  power  at  the 
wheels  (watts) 
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In  summary  the  variable  PW  is  calculated  as  a function 
of  W,  ACC,  HW  and  GRADE.  Subroutine  VEHDYN  does  this 
computation . 

3.2  VEHICLE  MOTOR 

In  this  section,  the  real  power  at  the  wheels  is  used 
to  calculate  the  voltage  and  current  at  the  motor  terminals. 
These  calculations  model  the  motor  and  drive  train  components 
in  the  propulsion  system. 

The  motor  is  a separately-excited  dc  machine  rated  at 
120  horsepower  and  600  volts  dc.  It  is  connected  to  the 
wheels  via  a gearbox  and  differential  whose  efficiency  is 
assumed  constant. 

Four  different  modes  of  operation  are  considered: 

1.  Motoring  with  voltage  control, 

2.  Motoring  with  field-weakening, 

3.  Regenerative  braking, 

4.  Friction  braking 

Modes  2 and  3 are  used  at  higher  speeds  - those  above 
a 30  mph  threshold.  Modes  1 and  2 are  used  if  the  power  at 


the  wheels  is  positive. 

More  precisely 

the  table  below  shows 

the  operating  conditions 

for  each  mode 

(W  is  velocity  (mph)  , 

PW  is  power  at  the  wheels 

(watts ) ) . 

MODE 

VELOCITY 

POWER 

1 

W<30 

PW>0 

2 

W>30 

PW>0 

3 

VV>30 

PW<0 

4 

VV<30 

PW<0 
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The  power  output  at  the  dc  machine  shaft  is  a function  of 
the  drive  train  efficiency  and  the  power  at  the  wheels. 

PW  = power  at  wheels  (watts) 

GREFF  = 0.92  gearbox  efficiency  (n.d.) 

PMS  = PW/GBEFF  power  at  machine  shaft  (watts)  if  mode  1 or  2 
PMS  = PW*GBEFF  power  at  machine  shaft  (watts)  if  mode  3 or  4 
Power  losses  in  the  dc  machine  are  either  stray  load  losses 
or  windage  losses. 

PMS  = power  at  machine  shaft  (watts) 

PL3  = 0.01  ABS(PMS)  where  PL3  are  the  stray  load  loss 

(watts)  computed  as  one  percent  of  the  absolute  value 
of  the  power  at  the  machine  shaft. 

W = vehicle  velocity  (mph) 

PL4  = 1461. 2*(W/30. ) **2.  5 where  PL4  are  the  windage 
losses  (watts) 

Note  PL4  is  modeled  as  a single  term  with  the  vehicle  velocity 
raised  to  the  2.5  power.  For  mode  1,  the  terminal  voltage  and 
current  are  calculated  as  follows: 

CF  = 10  where  CF  is  the  rated  field  current  (amp) 

C0N3  = 1.92  where  CON3  is  a motor  constant  (volts  per 
mph  per  amp) 

VV  = vehicle  velocity  (mph) 

EA  = C0N3*V*CF  where  EA  is  the  armature  back  EMF  (volt) 

PL3  = stray  load  losses  (watt) 

PL4  = windage  losses  (watt) 

VBD  = 2.0  the  comutator  brush  drop  (volts) 
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PMS  = power  at  motor  shaft  (watts) 

EA  = armature  back  EMF  (volts) 

CA  = (PMS  + PL3  + PL4)  / (EA  - VBD) 

The  terminal  voltage  can  then  be  found 
RA  = 0.15  armature  resistance  (ohms) 

VT  = EA+RA*CA  where  VT  is  the  terminal  voltage  (volts) 

For  MODE  2,  which  is  field  weakening,  the  terminal  voltage 
is  held  relatively  constant  at  its  30  mph  value,  while  the  field 
current  is  reduced  for  speed  control  above  30. 

Writing  first  an  equation  for  power-in  equals  power-out 
plus  losses , 

CA  = armature  current 

VT  = 576  rated  terminal  voltage  (volts) 

VBD  = 2.0  brush  drop  (volts) 

PMS  = power  at  machine  shaft  (watts) 

PL3  = strong  load  loss  (watts) 

PL4  = windage  loss  (watts) 

RA  = armature  resistance  (ohms) 

CA*VT  = PMS  + CA*  VBD  + CA**  2 * RA  + PL3  + PL4 
then  we  can  solve  it  for  the  armature  current. 

CA  = ((VT  - VBD)  + SQRT  ( (VT  - VBD)**  2+4*  RA*  (PMS+PL3+PL4) 
/ (2*RA) 

The  MODE  3,  which  is  regenerative  braking,  the  dc  machine 
acts  as  a generator.  In  this  case  the  armature  voltage  is 
greater  than  the  machine  terminal  voltage,  so  that  reverse 
armature  current  flows.  Figure  3-1  shows  this. 
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FIGURE  3-1.  THE  DC  MACHINE'S  ARMATURE  EQUIVALENT  CIRCUIT 
DURING  MOTORING  AND  REGENERATION. 
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In  this  case  the  armature  voltage  and  current  are  calculated 
as  shown: 

CF  = 10  rated  field  current  Camp) 

CON3  = 1.92  motor  constant  (volts  per  mph  per  amp) 

W = vehicle  velocity  (mph) 

EA  = CON3*W*CF  armature  back  EMF  (volts) 

PMS  = power  output  at  machine  shaft  (watts)  - note  PMS=0. 

PL3  = stray  load  loss  (watts) 

PL4  = windage  loss  (watts) 

CA  = (PMS  + PL3  + PL4)  / (EA  + VBD) 

RA  = Armature  resistance  (ohms) 

VT  = EA  + CA*RA 

For  MODE  4,  which  is  friction  braking,  the  motor  is  dis- 
connected and  mechanical  brakes  dissipate  the  energy.  The 
voltage  .and  current  are  set  essentially  to  zero. 

VT  = 0.01  machine  terminal  voltage  (volts) 

CA  = 0.01  machine  terminal  current  (amps) 

These  computations  are  carried  out  by  subroutine  MOTOR. 

The  final  step  in  MOTOR  is  to  calculate  the  magnitude  of  VT 
and  CA. 

VMTR  = ABS  (VT) 

IMTR  = ABS  (CA) 

TMTR  = 0 IF  CA>0 
TMTR  = 180  IF  CA<0 

In  summary  the  function  of  the  motor  simulation  is  to 
calculate  VMTR,  IMTR,  and  TMTR  using  PW  (the  power  at  the  wheels) 
and  W (the  vehicle  velocity) . 
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3.3  VEHICLE  POWER  - CONDITIONING  UNIT 


In  this  section,  the  equations  used  to  model  the  power- 
conditioning unit  (PCU)  are  discussed.  The  PCU  is  a 3-phase 
fully- controlled  bridge  rectifier,  which  interfaces  the  3- 
phase  power  rail  with  the  dc  machine  terminals.  When  real 
power  flows  from  the  power  rail  to  the  machine  terminals 
(Mode  1 or  2),  the  PCU  provides  a positive  terminal  voltage 
as  current  flows  against  the  machine's  back-EMF. 

When  real  power  flows  from  the  machine  terminals  to  the 
power  rail  (Mode  3),  the  machine's  back-EMF  is  in  the  opposite 
direction  from  above.  The  machine's  field  current  is  reversed 
to  produce  this  back  EMF  change.  The  machine's  back-EMF  is 
larger  than  the  PCU  output  (which  changes  size  as  the  firing 
angle  is  increased),  so  the  armature  current  flows  in  the 
same  direction.  Figure  3-3  shows  the  principle  components 
of  the  rectifier.  It  is  a six-pulse,  bridge  circuit  with  an 
isolation  transformer. 

The  analysis  here  follows  one  used  by  Schaefer  (Schaefer, 
1965).  The  SCR's  in  the  rectifier  have  a commutation  time  in 
which  the  current  flowing  through  one  SCR  is  switched  to 
another.  A voltage  drop,  EX,  is  associated  with  this  commuta- 
tion. It  is  a function  of  the  leakage  reactance  (on  the 
secondary)  of  the  transformer  and  the  output  current. 

CA  = output  current  (amps) 

XC  = 0.057  transformer  leakage  reactance  on  secondary  (ohms) 

EX  = 3 XC  *CA/3.1416  voltage  drop  due  to  commutation  (volts) 
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FIGURE  3-2.  THE  PCU  AND  DC  MACHINE  OPERATION 
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The  phase  delay  angle  or  thyristor  firing  angle  controls 
the  output  voltage.  The  following  equation  calculates  this 
angle,  ALPHA,  as  a function  of  the  commutation  voltage  drop, 
the  desired  output  voltage,  and  the  uncontrolled  output 
voltage . 

EX  = commutation  voltage  drop 

VRAIL  = complex  phase  voltage  at  the  power  rail  (p.u.) 
VBASE  = base  phase  voltage  (volts) 

EDO  = 3 SQRT  (6.)/3.1416  CABS  (VRAIL)  VBASE  the 
uncontrolled  output  dc  voltage  (volts) 

Note  SQRT(X)  is  square  root  of  X,  CABS(X)  is  the 
absolute  value  of  the  complex  quantity  X 
ED  = VT  the  desired  output  dc  voltage 
ALPHA  = ACOS  ((EX  + ED) /EDO)  the  SCR  firing  angle 
(radians ) 

Note  ACOS(X)  is  arc  cosine  of  X 

Using  this  firing  angle,  we  proceed  to  calculate  the 
power  into  the  PCU.  The  calculation  involves  a correction 
factor  to  account  for  the  effects  of  commutation.  For  this 
factor,  the  commutation  voltage  drop  is  converted  into  a 
commutation  angle. 

ED  = PCU  output  dc  voltage  (volts) 

EX  = commutation  voltage  drop  (volts) 

EDO  = uncontrolled  output  voltage  (volts) 

ALPHA  = firing  angle  (radians) 

AMU  = ACOS ( (ED-EX)/EDO)  - ALPHA  the  commutation  angle 
(radians ) 
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Note  ACOS (X)  is  arc  cosine  of  X 
The  correction  factor  is  calculated  in  three  steps: 

ALPHA  = firing  angle  (radians) 

AMU  = commutation  angle  (radians) 

CFl  = . 5 * (COS (ALPHA)  + COS (ALPHA  + AMU)) 

CF2  = (2.*AMU  + SIN(2*ALPHA)  - SIN(2*ALPHA  + AMU))/ 

(2. * (COS (ALPHA)  - COS (ALPHA  + AMU)) 

CFF  = 1,/SQRT(CF1**2  + CF2**2))  the  final  correction 
factor  for  commutation 
The  power  factor  is  next  calculated. 

ALPHA  = firing  angle  (radians) 

CFF  = commutation  correction  factor 

PF  = COS (ALPHA) *CFF  the  power  factor  at  the  PCU  input 
The  real  power  into  the  PCU  is  equal  to  the  real  power, 
assuming  negligible  losses. 

VT  = PCU  output  voltage  (volts) 

CA  = PCU  output  current  (amps) 

REF  = VT*CA  real  power  input  to  the  PCU  (watts)  in  all 
three  phases. 

The  reactive  power  into  the  PCU  is  calculated  below. 

EDO  = uncontrolled  output  voltage  (volts) 

CF2  = correction  failure  for  ratio  of  fundamental 
reactive  power  to  uncontrolled  dc  power 
CA  = output  current 

AIMP  = ED0*CA*CF2  the  reactive  power  (VA)  in  all  three 
phases . 
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The  final  equations  are  concerned  with  adding  in  the 
auxiliary  loads  and  then  normalizing  the  power  values  to 
per  unit. 

PAUX  = 11,520  the  real  power  for  the  auxiliary  load 
(watts)  in  all  3 phases. 

QAUX  = 8,640  the  reactive  power  for  the  auxiliary 
loads  (VA)  in  all  3 phases. 

PEASE  = 333,000  BASE  POWER  (VA)  for  one  phase 
BP  = (REF  + PAUX) /PBASE/ 3 per  unit  real  power  at  the 
power  rail/vehicle  interface. 

BQ  = (AIMP  + QAUX) /PBASE/ 3 . per  unit  reactive  power  at 
the  power  rail/vehicle  interface. 

3.4  ADMITTANCE  MATRIX 

The  power  rail  segments  and  power  cables  each  have  an 
impedance  which  affects  the  power  flow  from  the  utility 
connection  to  the  vehicles.  This  impedance  is  modeled  as  a 
linear  function  of  the  guideway  segment  length  or  as  a 
constant  for  the  power  cables  (and  step-down  transformers). 
The  impedance  values  are  calculated  as  per  unit  impedances, 
converted  to  per  unit  admittances,  and  stored  in  a sparse 
admittance  matrix.  The  admittance  computation  is  summarized 
below . 

CASE  1.  A POWER  RAIL  BETWEEN  NODE  I AND  J 

ZGY  = complex  power  rail  impedance  (per  unit  ohms  per  ft 
per  phase) 
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D = distance  along  the  power  rail  (ft)  from  node  I to 
node  J 

YD  = l./(D*ZGY) 

YD  = complex  power  rail  admittance  (per  unit  mhos 
per  phase)  between  nodes  I and  J 
CASE  2.  A CABLE  OR  TRANSFORMER  BETWEEN  NODE  I AND  J. 

GWAY(II,4)  = real  part  of  the  complex  fixed  impedance. 

between  nodes  I and  J (see  Section  6.  using  the  model.) 
GWAY(II,5)  = imaginary  part  of  the  complex  fixed  impedance 

between  nodes  I and  J.  (See  Section  6.  using  the  model.) 
YD  = l./CMPLX(GWAY(II,4) , GWAY(II,5) 

YD  = complex  power  rail  admittance  (per  unit  mhos  per 
phase)  between  nodes  I and  J. 

The  admittance  computation  involves  nodes  in  the  guideway 
network.  A numbering  convention  is  adopted  for  these  electrical 
nodes.  Node  1 is  the  electric  utility  connection.  Nodes  2, 

3,  ...,  NCAR+1  are  vehicles. 

In  summary,  the  function  of  the  PCU  simulation  is  to 
calculate : 

1)  BP  - per  unit  real  power  into  PCU  at  power  rail/ 

vehicle  interface 

2)  BQ  - per  unit  reactive  power  into  PCU  at  the  power 

rail/vehicle  interface 

3)  ALPHA  - SCR  firing  angle  (radians),  given  as  inputs. 
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1)  UMTR  - magnitude  o£  the  motor  terminal  voltage  (volts) 

2)  IMTR  - magnitude  of  the  motor  terminal  current  (amps) 

3)  TMTR  - has  value  0 if  motor  terminal  voltage  is 

greater  than  0,  has  value  180  otherwise 

4)  VRAIL  - complex  power  rail  voltage  (volts) 

This  simulation  is  carried  out  by  the  subroutine  PCU. 

These  nodes  are  usually  moving  and  have  a current  associated 
position  (a  guideway  segment  and  displacement  from  the  segment 
beginning).  Nodes  NCAR+2,  NCAR+3 , . . . ,SZ-1-1,SZ1  are  fixed 
nodes  along  the  guideway.  These  nodes  include  power  feed 
points,  diverge  or  merge  power  rail  connections  and  power  tie 
points  between  parallel  segments.  There  are  SZl  nodes  in 
total  where  SZl  is  a user- supp lied  parameter. 

The  admittance  matrix  is  stored  in  several  tables,  but 
only  nonzero  values  are  stored  in  order  to  save  space.  The 
tables  are  shown  in  Figure  3-4  and  their  entries  are  described 
below.  Note  that  in  the  two-dimensional  admittance  row  I and 
column  I are  devoted  to  node  number  I.  The  entry  Y(I,I)  is 
the  driving  point  admittance  at  node  I.  The  entry  Y(I,J), 
where  I f J,is  the  transfer  admittance.  Finally  since  Y(I,J) 
equals  Y(J,I),  the  tables  of  Figure  3-4  have  only  one  entry 
for  those  pairs.  Again,  this  is  done  to  save  space. 

3.5  LOAD  FLOW  ALGORITHM 

The  equivalent  circuit  for  the  power  distribution  network 
is  described  by  a set  of  simultaneous  equations.  They  relate 
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FIGURE  3-4.  a)  AN  ADMITTANCE  MATRIX  AND,  b)  ITS  TABULAR 
REPRESENTATION 
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the  voltages,  current,  and  power  at  each  node  in  the  network. 

The  load  flow  algorithm  is  just  the  iterative  solution  tech- 
nique used  on  the  equations.  Many  variations  exist  (Stott, 
1974),  but  the  one  used  here  is  a power-mismatch  version  in 
which  a generalized  Newton-Raphson  method  is  used.  The  basic 
equations  are  given  below.  All  quantities  are  fundamental 
(60HZ)  values . 

DP(k)  - jDQ(k)  = CONJG(E(k))*I (k) 

DP(k)  = real  power  (per  unit)  delivered  into  the  network 
at  node  k 

DQ(k)  = reactive  power  (per  unit)  delivered  into  the  net- 
work at  node  k 

E(k)  = line  to  neutral  voltage  (per  unit)  at  node  k 

I (k)  = phase  current  (per  unit)  into  the  network  at  node  k. 

CONJG(E(k))  = complex  conjugate  of  E(k) 

The  admittance  values  are  brought  into  the  analysis  with  the 
following  network  performance  equation: 

SZl 

I (k)  = Z!  y (k,n) *E(n) 
n=  1 

I (k)  = phase  current  (per  unit)  into  the  network  at 
node  k 

y(k,n)  = admittance  (per  unit)  value  from  the  admittance 
matrix.  It  is  either  a driving  point  admittance 
(if  k=n)  or  a transfer  admittance  (if  k^n) . 

E(n)  = line  to  neutral  voltage  (per  unit)  at  node  n. 

SZl  = number  of  network  nodes 
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By  combining  these  two  equations  the  following  set  of 
nonlinear  equations  results: 

SZl 

DP(k)  - jDQ(k)  = CONJG(E(k))  X)  y(k,n)*E(n) 

n=l 

The  solution  involves  the  following  steps: 

1.  Find  the  admittance  matrix  so  that  y(k,n)  values  are 
known. 

2.  Assume  initial  voltages  (E(k))  at  each  node. 

3.  Find  the  P(k)  and  QCk)  values  at  each  node  in  the 
network  where  a vehicle  needs  electric  power.  Note 
that  the  real  power  could  be  positive  to  model  a 
regenerating  vehicle  or  negative  to  model  a power 
absorbing  vehicle.  These  values  are  the  scheduled  power. 
They  are  functions  of  the  vehicle's  propulsion  system, 
its  current  mission  profile  demand,  and  the  node 
voltage  E(n)  at  the  power  rail/vehicle  interface. 

4.  Using  the  nonlinear  system  of  network  equations, 
with  the  admittance  values  y(k,n)  and  node  voltages 
E(n),  find  the  delivered  power.  Any  mismatch  between 
the  scheduled  and  delivered  power  can  be  corrected  by 
revising  the  voltages  ECn)  at  each  node.  Succeeding 
steps  perform  this  correction. 

5.  Calculate  the  differences  between  scheduled  and 
delivered  power. 
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DELTAP  (k)  - P(k)  - DP(k) 

DELTAQ  (k)  = Q(k)  - DQ(k) 

P(k)  = real  power  scheduled  (per  unit)  into  the 
network  at  node  k 

DP(k)  = real  power  delivered  (per  unit)  into  the 
network  at  node  k 

Q(k)  = reactive  power  delivered  (per  unit)  into  the 
network  at  node  k 

DQ(k)  = reactive  power  delivered  (per  unit)  into  the 
network  at  node  k 

DELTAP (k)  = real  power  (per  unit)  mismatch  at  node  k 
DELTAQ(k)  = reactive  power  (per  unit)  mismatch  at 
node  k 

6.  The  iteration  may  stop  if  the  maximum  DELTAP (k)  and 
DELTAQ(k)  are  below  some  cutoff,  EPSIL  (Typically  0.10) 
If  the  maximum  mismatch  is  too  large,  then  the 
following  steps  revise  the  node  voltages. 

7.  Calculate  the  current  into  the  network  at  each  node 

I(k). 

I(k)  = (P(k)  - jQ(k))/CONJG(E(k)) 

8.  Solve  a set  of  linear  equations  for  the  change  in 
node  voltages.  These  equations  are  shown  in  Figure 
3-5.  The  Jacobian  matrix,  J,  is  a sparse  matrix,  so 
only  the  nonzero  elements  are  stored  in  a table. 

Figure  3-6  depicts  these  tables.  The  standard  Gauss- 
Jordan  reduction  algorithm  is  used  to  find  the  updated 
E(k)  values. 
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FIGURE  3-7.  EQUIVALENT  CIRCUIT  FOR  ONE  PHASE  OF  A POWER  RAIL 
SEGMENT.  THE  IMPEDENCE  Z IS  A FUNCTION  OF  LENGTH 
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9.  The  computation  continues  at  step  3. 

This  computation  is  carried  out  by  the  subroutine  named 
NR.  In  addition  the  subroutine  FORMY  is  used  for  the 
admittance  matrix  formation  and  subroutines  GETJ  and  PUTJ 
handle  the  Jacobian  matrix  manipulation. 


3.6'  DISTORTION  POWER 

In  this  section  the  high  frequency  currents  are  used 
to  calculate  the  power  components  which  are  not  included  in 
the  fundamental  real  and  reactive  power.  The  analysis  of  the 
high  frequency  currents  is  concerned  with  two  types  of  loca- 
tions in  the  power  distribution  network.  One  is  the  power  rail/ 
vehicle  interface  where  the  nonlinear  PCU  generates  the 
harmonic  currents.  Second  is  the  utility  connection  where  the 
currents  impact  the  supply  network  and  where  the  metering  is 
located . 

In  order  to  simplify  the  analysis,  two  assumptions  are 
made:  (1)  the  voltage  waveform  at  the  utility  connection  is 

sinusoidal  (i.e.,  it  is  a stiff  source),  and  (2)  the  power 
rail  is  a short  transmission  line  (i.e.,  it  is  a constant 
series  impedance  dependent  on  rail  length) . The  equivalent 
circuit  for  a power  rail  segment  is  shown  in  Figure  3-7. 

This  analysis  involves  two  steps.  First  the  harmonic 
currents  at  each  power  rail/vehicle  interface  are  calculated. 
Second  these  currents  are  combined  into  the  current  at  the 
utility  connection. 
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FIGURE  3-8.  LINE  TO  NEUTRAL  VOLTAGE  AT  THE  UTILITY  AND  AT 
POWER  RAIL/VEHICLE  INTERFACE 


FIGURE  3-9-  LINE  CURRENT  AT  THE  POWER  RAIL/VEHICLE  INTERFACE 
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The  harmonic  currents  are  calculated  by  Fourier  analysis 
on  the  wave  forms.  The  current  wave  form  used  is  shown  in 
Figure  3-8.  This  wave  form  is  discussed  in  the  literature 
(Schaefer,  1965,  p 327-334).  It  represents  the  line  current 
for  one  phase  when  two  effects  are  present:  (1)  a delay  due 
to  the  SCR  firing  angle  a,  and  (2)  a rise  (fall)  which  is  not 
instantaneous  commutation. 

The  expression  to -describe  the  current  curve  is  given 
below  for  the  three  segments  (rise,  constant,  fall), 
a = SCR  Firing  Angle  (radians) 

Y = 2it/3  Length  of  current  pulse  with  instantaneous 
commutation  (radians) 

U = commutation  angle  (radians) 

ID  = output  dc  current  from  PCU  (amp) 

IB  = per  unit  base  Current  (amp) 

IR  = PCU  input  line  current  during  Rise  Time  as  a function 
of  cot  (per  unit  amp) 

IR(cot)  = ID*(cos(a)  - cos  (cot+ y/ 2)  ) / (cos  (a)  - cos  (a+y)  ) /IB 
IC(cot)  = ID/IB  PCU  input  line  current  during  constant  time 
as  a function  of  cot  (per  unit  amp) 

IF  = PCU  input  line  current  during  fall  time  as  a function 
of  cot  (per  unit  amp) 

IF(cot)  = I D*  (1- (cos  (a)  - cos  (cot- y/ 2)  ) / (cos  (a)  - cos  (a+y) )) /IB 
The  Fourier  analysis,  which  yields  the  harmonic  currents,  is 
performed  next.  Because  of  the  positive/negative  symmetry  in 
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the  current  wave  form,  the  integral  is  evaluated  over  a half 
period. 

ID  = output  dc  current  from  PCU(amp) 
a = SCR  firing  angle  (radians) 
y = commutation  angle  (radians) 

IR((jot)  = PCU  input  line  current  during  rise  time  (per  unit  amp) 
IC(cjt)  = PCU  input  line  current  during  constant  time  (per 
unit  amp) 

IF(o3t)  = PCU  input  line  current  during  fall  time  (per  unit  amp) 
A(N)  = (SQRT(2))/tt* 


/ 


a+y-y/2 

IR(o)t)  *sin  (Nwt)  dtut 


a-y/2 


a+y/2 

^ I C (cot)  *sin  (Nwt)  dtot 
a+y-y/2 

a+y-y/2 

J'  IF(o)t)  *sin  (Nwt)  dwt 
a+y/  2 

A(N)  = RMS  value  of  the  N-th  line  current  harmonic 
the  quadrature  component  (per  unit  amp). 


B(N)  = (SQRT(2)  )/tt* 


a+y-y/2 

^ IR(a)t)*cos(Noot)  dcot 
a-y/2 

a+y/2 

J'  IC(oat)*cos  (No)t)  dcot 
a+y-y/2 

a+y+y/ 2 

J'  IF(cot)*cos(Ncot)  dcot 
a+y/2 


(B/N)  = RMS  value  of  the  N-th  line  current  harmonic-in  phase 
component  (per  unit  AMP).  Appendix  C shows  the  result  of 
evaluating  these  integrals. 

The  final  step  in  the  analysis  is  the  combining  of  the 
harmonic  current  components  at  the  power  rail/vehicle  inter- 
face into  harmonic  currents  at  the  utility  connection.  For 
this,  each  order  harmionic  is  considered  separately.  So 
in  the  discussion  that  follows  the  fifth  harmonic  (N=5)  is 
used  but  N could  be  5 , 7 , 11 , 13 , . . . , 6N- 1 , 6N+1 . Figure  3-9 

shows  a small  power  rail  distribution  system. 

t h 

By  using  superposition,  the  harmonic  current  at  the 

utility  connection,  1^^,  can  be  found.  Each  current  source 
corresponding  to  a vehicle  is  treated  separately  as  the  others 
are  replaced  by  their  infinite  internal  impedance.  Since  , 
Z2,  and  Z^  are  series  impedances  for  the  short  transmission 
line  model,  the  current  is  I2  + I-^.  In  general  the 
harmonic  current  flowing  from  the  utility  into  the  power  rail 
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network  is  the  sum  of  the  harmonic  currents  flowing  into 

the  power  rail/vehicle  nodes.  The  sum  is  a vector  sum  taking 
into  account  the  phase  angles.  To  illustrate,  Figure  3-10 
is  presented. 

Each  fifth  harmonic  at  a vehicle  has  a phase  angle  shown 

in  Figure  3-10  as  <p , which  is  the  number  of  degrees  by  which  it 

lags  the  fundamental  voltage  at  the  utility  (our  reference). 

The  angle  tj)  is  calculated  by  the  following  equation: 

N = harmonic  number  (n.d.) 

th 

B(N)  = RMS  value  of  the  N line  current  harmonic-in 

phase  component  (per  unit  amp) . 

t h 

A(N)  = RMS  value  of  the  line  current  harmonic-quadrature 

component  ( per  unit  amp) 

(p2  = ATAN(A(N)/B(N) 

th 

(p2  = angle  by  which  N harmonic  current  at  vehicle  lags 
fundamental  voltage  at  vehicle  (radians) 

(pi  = angle  by  which  fundamental  voltage  at  the  vehicle 
lags  fundamental  voltage  at  the  utility  (radians) 
cj)  = N * (pi  + (p2 

t h 

(})  = angle  by  which  the  N harmonic  current  at  the 

power  rail/vehicle  interface  lags  the  fundamental 
voltage  at  the  utility  (radians) 

The  fifth  harmonic  at  the  electric  utility  is  only  one  of 
several  harmonics  which  may  be  present.  The  RMS  values  of  all 
these  harmonics  can  be  found  from  the  sum  of  the  squares. 
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CURE  3-11.  HARMONIC  CURRENT  AT  THE  POWER  RAIL/ VEHICLES 
WITH  VOLTAGE  WAVEFORMS  FOR  REFERENCE 
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= RMS  value  of  the  harmonic  current  at  the 

utility  node 

DC(1)  = SQRT  (E  I(N)**2  ) 

N 

DC(1)  = RMS  value  of  the  distortion  current  at  the 
utility  node  (per  unit  amp) . 

The  distortion  power  can  be  calculated  as  the  product  of  the 
utility  voltage  and  distortion  current. 

DC(1)  = RMS  distortion  current  (p.u.) 

E(l)  = utility  voltage  (L-N)  (p.u.) 

DKVA  = E(l)  * DC(1) 

DKVA  = distortion  power  at  the  utility  (p.u.) 

The  total  apparent  KVA  may  also  be  calculated. 

P(l)  = active  power  at  the  utility  (p.u.) 

Q(l)  = reactive  power  at  the  utility  (p.u.) 

DKVA  = distortion  power  at  the  utility  (p.u.) 

TKVA  = SQRT(P(1)**2+Q(1)**2+DKVA**2) 

TKVA  = total  apparent  KVA  at  the  utility  (p.u.) 

The  distortion  power  is  calculated  from  a Fourier  analysis 
of  the  current  wave  shape.  The  parameters  of  the  waveshape  are 
a (SCR  firing  angle)  and  p (commutation  angle) . The  outputs  are 
the  distortion  power  and  total  apparent  KVA  at  the  electric 
utility  connection  node.  This  modeling  is  done  by  routine 
DISTOR. 

3.7  VEHICLE  POSITION  UPDATES 

Periodically  during  the  simulation,  the  vehicle  positions 
are  updated.  This  section  presents  the  analysis  used  to  calculate 
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these  updated  positions.  The  analysis  involves  data  for  each 
vehicle  from  three  sources: 

1.  The  commanded  mission  profile  point,  which  is  a 
requested  velocity,  and  encountered  grade,  encountered  head- 
wind, a position  on  the  guideway  and  a time  interval  before 
the  next  profile  point  should  be  attained. 

2.  The  current  vehicle  state,  which  is  like  the  commanded 
mission  profile  point  but  with  actual  rather  than  commanded 
values  used. 

3.  The  propulsion  system  limitations  for  the  maximum 
acceleration,  jerk,  thrust  and  power  at  the  wheels. 

The  strategy  for  using  this  data  is  illustrated  in  Figure  3-11. 
Three  commanded  mission  profile  velocity  points  (CV(t),  (CV(t+l) , 
CV(t+2))  are  shown,  together  with  the  actual  velocity  (AV(t+l)). 
Each  commanded  velocity  is  treated  as  a speed  limit  (which  the 
vehicle  accelerates  or  brakes  to  achieve) . These  limits 
change  in  two  ways.  First,  if  a vehicle  passes  the  guideway 
position  associated  with  a commanded  velocity  profile  point, 
the  commanded  velocity  changes.  That  is,  referring  to  Figure 
3-12,  each  point  in  time  has  a position  (P(t),  P(t+2)).  It 
is  assumed  that  the  guideway  conditions  (grade,  headwind,  speed 
limit)  which  are  recorded  at  a time  such  as  t+1  for  position 
P(t+1)  remain  constant  until  position  P(t+2).  A problem 
arises  when  a vehicle  is  dwelling  (velocity  is  zero) . In 
this  case  it  will  never  pass  the  next  guideway  position  and 
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FIGURE  3-12. 


Ca)  COMMAND  AND  ACTUAL  VELOCITY  PROFILE 
(b)  POSITION  PROFILE  POINTS. 


''move”  to  a new  commanded  mission  profile  point.  To  correct 
for  this  a second  method  is  employed  to  change  mission  profile 
points.  This  one  says  if  the  velocity  is  zero  and  the  time 
interval  is  exceeded,  then  use  the  next  profile  point.  The 
next  point  should  have  a nonzero  velocity,  which  is  used  to 
calculate  an  acceleration  to  move  the  vehicle.  The  position 
of  each  vehicle  is  calculated  by  the  following  formula: 

DT  = time  interval  (sec) 

POS  = vehicle  displacement  (ft)  from  last  guideway 

reference  point.  (Its  position  at  some  previous 
time  t.) 

WEL  = vehicle  velocity  (mph)  at  time  t (actual  velocity 
achieved) 

VACC  = vehicle  acceleration  (mph/sec)  during  interval  DT 
Cl  = 1.4667  conversion  factor  to  change  mph  into 
feet/sec . 

NEWPOS  = POS  + WEL  + DT  * Cl  + 0.5  *VACC  * Cl  * DT  ^ DT 
NEWPOS  = vehicle  displacement  (feet)  from  last  guideway 
reference  point.  (Its  position  at  time  t + DT) 

It  should  be  noted  that  the  acceleration  VACC  takes  into 
account  jerk  acceleration  and  thrust  limitations.  The  new 
velocity  can  also  be  found. 

DLEVEL  = vehicle  velocity  (mph)  at  time  t 

VACC  = vehicle  acceleration  (mph/sec)  during  interval  DT 
DT  = time  interval  (sec) 

WEL  = vehicle  velocity  (mph)  at  time  t + DT 
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This  simulation  is  done  by  subroutine  NUPNT  which  uses  the 
mission  profile  data,  current  vehicle  state  information,  and 
subroutine  VEHDYN  to  calculate  a new  position  and  velocity 
for  the  entire  vehicle  fleet. 


54 


4.  SAMPLE  OUTPUT  - AGRT 


In  this  section,  a sample  run  is  described.  From  this 
run  the  reader  can  see  the  type  of  output  data  produced  by 
the  simulation.  This  data,  together  with  the  mission  pro- 
file input  data  (Appendix  A) , provides  a complete  check  case 
for  the  program. 

The  output  consists  of  three  parts:  (1)  a record  of  the 

guideway  data  used  for  input,  (2)  a record  of  the  mission 
profile  data  used  for  input  (it  is  a partial  listing  of  only 
the  first  15  lines  of  the  table) , (3)  output  tables 
showing  the  vehicles’  positions,  the  power  flow  from  the 
utility  connection,  the  network  node  voltages  , and  the  state 
of  the  propulsion  system  for  each  vehicle.  The  latter  group 
(3)  of  output  tables  are  repeated  for  each  sample  time. 

Table  4-1  is  the  complete  output  for  two  vehicles  run- 
ning along  a loop  (similar  to  the  Morgantown  system) . The 
power  rail  is  pictured  in  Figure  4-1.  In  this  case  the  two 
vehicles  become  electrical  nodes  (buses)  numbered  2 and  3. 

The  reader  should  note  that  the  vehicles  are  identified  as 
vehicle  number  1 , 2 , . . . , NCAR,  but  once  they  are  included  in 
the  equivalent  electric  circuit,  they  are  node  (bus)  number 
2,3,... ,NCAR+1. 

The  list  below  describes  the  output  data: 

1.  Guideway  segments  at  time -n  sec.  This  table  dis- 
plays the  guideway  input  data  but  appends  the  node 
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number  for  any  vehicle  on  that  segment. 

2.  Power  flow  summary.  This  table  displays  the  power 
flow  along  each  individual  cable  which  is  connected 
to  the  electric  utility  connection.  It  is  funda- 
mental P and  Q in  per  unit. 

3.  Total  apparent  KVA.  This  is  the  apparent  power 
supplied  by  the  utility  in  per  unit. 

4.  Bus  voltages. 

This  table  lists. the  bus  voltage  in  per  unit,  the 
fundamental  power  into  the  network  in  per  unit 
(loads  are  negative,  generators  are  positive)  and 
the  distortion  current  in  per  unit. 

5.  State  of  the  vehicles 

For  each  vehicle  number  (1  to  NCAR)  this  table  lists 
a vehicle  position  (a  guideway  identifier  and  a 
footage  displacement) , current  velocity  (in  MPH) , 
acceleration  (in  MPH/s)  ideal  or  commanded  velocity 
(in  MPH) , encountered  grade  (in  percent) , encountered 
headwind  (in  MPH) , motor  terminal  voltage  (volts) , 
motor  terminal  current  (amps)  number  of  passengers, 
number  of  cars  coupled  together,  a mode  (1-voltage 
controlled  speed,  2 - field-weakening  controlled  speed, 
3-regenerative-braking,  4-friction  braking),  the 
power  conditioning  units  firing  angle  (alpha-in  ra- 
dians) . 

These  tables  are  written  by  the  subroutines  named:  RITE 

or  VEHSUM.  The  labels,  formats  and  variables  are  described  there. 
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TABLE  4-1. 


HORGANTOWN  GUIDEWAV  DATA  - THREE  LINES 


GUIOEWAV 
FROM  TO 

DATA 

SEGMENT  START 

END 

6 • 6* 

6 • S • 

E 11650.000 

M3  12295.000 

12285.000 
17  160.000 

S.  4. 
A. 

M4  17160.000 

W 1.000 

21 903.000 
2050.000 

4.  5. 

5.  6. 

Ml  2050.000 

M2  6785.000 

6785.000 

11650.000 

1.  4. 

1.  5. 

TT  0.000 

TT  0.000 

0.025 

0.025 

1 . 6. 

4GT  2 

TT  0.000 

0.025 

0 746 

0 385 

73  224  0 

435  584  0 

2 2 
1 413 

0 0 
3.0000 

0.0000  0.0000 

0. oooow 

903.11 

1.0000 

1.0000 

1.0200  0.0000 
2.0400  0.0000 

o.oooow 
0. oooow 

905.35 

906.44 

1.0000 

1.0000 

2.1760  0.0000 

2.3120  0.0000 

o.oooow 

o.oooow 

911.73 

915.23 

1.0000 
I. 0000 

2.4480  0.0000 

2.5840  0.0000 

o.oooow 

o.oooow 

918.92 

922.81 

1.0000 

1.0000 

2.7200  0.0000 

4.0800  0.0000 

o.oooow 
0. oooow 

927.79 

934.77 

1.0000 

1.0000 

5.4400  0.0000 

5.4400  0.0000 

o.oooow 

o.oooow 

942.75 

950.73 

1.0000 

1.0000 

5.4400  0.0000 

5.4400  0.0000 

o.oooow 
0. oooow 

958.71 
9 66.69 

1 .0000 
1.0000 

5.4400  0.0000 

5.4400  0.0000 

o.oooow 

o.oooow 

9 74  .67 
982.65 

GUIDEWAY 
6.  6.E 

SEGMENTS  AT  TIME  « 0 SEC 

11650.000  12265.000  3. 

0.  0. 

0.  0. 

6.  5. M3 
5.  4.M4 

12285.000  17160.000  0. 

17160.000  21903.000  0. 

0.  0. 

0.  0. 

0.  0. 

0.  0. 

4.  4.W 
4.  5.M1 

1.000  2050.000  2. 

2050.000  6765.000  0. 

0.  c. 
0.  0. 

0.  0. 

0.  0. 

5.  6. M2 
1.  4.TT 

6785.000  11650.000  0. 

0.000  0.025  0. 

0.  0. 

0.  0. 

0.  0. 

0.  0. 

1.  5.TT 
1.  6.TT 

0.000 

0.000 

0.025  0. 

0.025  0. 

0.  0. 
c.  c. 

0.  0. 

0.  0. 

POWER  FLOW  SUMMARY  AT  SIMULATED  CLOCK 

TIME  OF 

0 SEC  : 

POWER  FLOW  FROM  NCOE  1 TO 
POWER  FLOW  FROM  NODE  1 TO 

NODE  4 IS 

NODE  5 IS 

0.0108 
0. 0044 

-0.0086 

-0.0030 

POWER  FLOW  FROM  NObE  1 TO 

NODE  6 IS 

6.014A 

-0.0540 

THE  TOTAL  APPARENT  KVA  IS 

0.073040 

IN  PER  UNIT 

BUS 

1 

VOLTAGE 

1.0000  0.0000 

POWER  SCHEDULED 
0.0300  0.0655 

POWER  DELIVERED 
0.0300  0.0655 

DISTORTION  CURRENT 
0.0119 

2 

3 

0.9990  -0.0005 

0.9961  -0.0001 

-0.0115 

-0.0184 

-0.0087 

-0.0570 

-0.0115 

-0.0184 

-0  .0  08  7 
-0.0  570 

0.0000 
0 .0119 

4 

5 

0.9998  -0.0003 

0.9999  -0.0001 

J.OOOO 

0.0000 

0.0000 

0.0000 

0.0000 

-0.0000 

0.0000 

0.0000 

0.0000 

O.OOCO 

I 
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TABLE  4-1  (CONI' D) 
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5.  VALIDATION  STUDIES 


In  order  to  validate  the  model,  two  tests  were  performed. 
The  first  compared  the  program's  output  to  other  simulations, 
while  the  second  compared  it  to  an  operating  system. 

For  the  first  test,  we  used  a model  of  the  Dallas-Fort 
Worth  Airport  AIRTRANS  system  developed  by  Vought  Corporation. 
Several  simulation  studies  had  been  done  using  this  model. 

In  essence,  the  AIRTRANS  system  consisted  of  fifteen  (one  or 
two  car)  trains  running  over  fourteen  miles  of  guideway.  Four 
check  cases  were  chosen,  and  the  results  of  the  Vought  simula- 
tion and  our  simulation  were  compared  for  each  case.  The  re- 
sults show  the  electric  power  values  to  be  similar.  Any 
differences  are  accounted  for  by  truncation  errors 
(Typically  II) . 

For  the  second  test,  the  four  check  cases  were  compared 
with  values  measured  at  the  Dallas  Fort-Worth  Airport.  Our 
check  cases  are  within  10  percent  of  the  typical  average  real 
power . 
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6.  USING  THE  MODEL 


Three  steps  are  involved  in  using  the  model:  (1)  set 

up  the  propulsion  system  subroutines,  (2)  create  a data  file 
for  the  guideway,  and  (3)  create  a data  file  for  the  mission 
profile . 

Step  one  requires  the  user  to  modify  the  analysis  pre- 
sented in  section  3 and  then  change  the  subroutines  which  em- 
bed that  analysis  in  the  simulation. 

Steps  two  and  three  involve  input  data  preparation. 

Two  data  files  are  used  by  the  model.  The  first  contains 
guideway  information;  the  second  contains  mission  profile  in- 
formation. The  size  of  these  data  files  is  recorded  in  com- 
pile-time parameters. 

First  - consider  the  guideway  file.  It  is  a series  of 
records  which  are  read  sequentially  into  the  array  GWAY.  The 
parameter  GWAYRO  specifies  the  number  of  records  or  lines. 

Each  record  is  read  into  a row  of  the  array  GWAY  at  columns 
1 to  5.  The  record  describes  one  segment  of  guideway  which 
connects  two  points  or  nodes  on  the  guideway.  The  following 
naming  convention  is  used  for  the  guideway:  Guideway  nodes 

are  integers.  Node  one  is  the  utility  connection.  Nodes 
two  and  higher  are  the  power  cross-unders  and  power  feed 
points  in  any  order. 
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The  following  list  describes  the  record  fields; 

1.  FN  - the  node  number  at  the  beginning  of  the 
segment.  Vehicles  enter  the  segment  at  this 
node . 

2.  TN  - the  node  number  at  the  end  of  the  segment. 
Vehicles  leave  the  segment  at  this  node. 

3.  SI  - the  two -character  segment  identifier.  This 
is  a unique  name  for  each  segment. 

4.  SD  - a distance  mark  for  the  FN  node.  Other 
points  along  the  guideway  are  measured  with 
respect  to  this  reference.  The  units  are  in 
feet . 

5.  ED  - a distance  mark  for  the  TN  node.  The  length 
(feet)  of  the  segment  is  ED-SD. 

By  convention  ED>SD. 

The  subrouting  SETGW  reads  in  these  records  and  stores 
them  in  array  GWAY.  The  compile -time  parameter,  GWAYRO,  de- 
fines the  number  of  rows  in  array  GWAY  so  that  the  number  of 
rows  equals  the  number  of  records. 

The  number  of  columns  in  the  array  GWAY  is  usually  10. 
The  first  five  columns  are  used  for  the  guideway  records  and 
the  remaining  are  used  for  storing  the  node  number  of  any 
vehicle  which  is  currently  positioned  on  that  guideway  seg- 
ment . 

One  exception  for  the  guideway  data  file  occurs  for  the 
case  in  which  a step-down  transformer  and/or  cables  (rather 
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FN 

TN 

SI 

SD 

ED 

2X 

3.0 

F3.0 

A2 

F6.0 

F6.0 

FIGURE  6-1  GUIDEWAY  RECORD  AND  FORMAT 

than  a power  rail)  connect  two  nodes.  In  that  case,  the  SI 
field  has  a special  value,  ”TT''.  The  transformer  impedance 
(R+jX  in  per  unit)  is  stored  in  the  SD  and  ED  fields  for  R 
and  X respectively. 

The  second  data  file  contains  mission  profile  information. 
It  is  a series  of  records  which  are  read  sequentially  into  the 
array  PFILE.  The  parameter  PFILER  specifies  the  number  of 
records.  Each  record  is  read  into  a row  of  the  array.  The 
record  describes  one  point  of  the  mission  profile.  An  arbitrary 
number  of  records  are  collected  together  and  called  a route 
segment.  Figure  6-2  shows  this.  The  data  at  each  mission 
profile  point  is  described  in  the  following  list: 
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A ROUTE  SEGMENT  IN  ARRAY  PFILE 


DELTA 

TIME  VEL  GRADE 

(SEC)  (MPH)  (%) 


1.0 

1.0 

1.0 

1.0 

1.0 

0.0 

4F9.4 


A2  F9.2 


FIRST  RECORD 


LAST,  RECORD 
FORMAT 


FIGURE  6-2  CONTIGUOUS  ROWS  IN  ARRAY  PFILE 
WHICH  FORII  A ROUTE  SEGMENT. 
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1.  DELTA  TIME  - A NORMAL  TIME  interval  (seconds,  be- 

tween this  point  and  the  succeeding 
profile  point.  The  last  record  in  a 
route  segment  has  DELTA  TIME  as  0. 

2.  VELOCITY  - The  COMMANDED  SPEED  UNIT  (MPH) 

3.  GRADE  - The  encountered  grade  (I) 

4.  HEADWIND  - The  encountered  headwind  (MPH) 

5.  LOCATION  CODE  - A segment  identifier  (two  unique 

characters)  and  a distance  (feet). 

Each  such  route  segment  is  identified  by  the  row  number  in 
PFILE  where  the  first  record  is.  Figure  6-3  shows  a collection 
of  route  segments  in  the  array  PFILE.  A complete  mission  pro- 
file data  file  is  contained  in  Appendix  A . 

A vehicle's  mission  is  usually  a sequence  of  route  seg- 
ments as  Figure  6-4a  indicates.  So  this  sequence  of  route 
segments  is  part  of  a prefix  which  is  included  in  the  mission 
profile  data.  The  following  list  describes  those  data  items: 

1.  Title  - a forty  character  heading  which  appears  on 
the  output 

2.  ICAR  - an  array  with  one  column  for  each  vehicle, 
and  an  arbitrary  number  of  rows.  Column  j contains 
the  route  segments  (by  first  row  number)  in  the 
mission  profile  of  vehicle  j . The  first  and  last 
segments  are  arbitrarily  zero. 

3.  SROWl-  a vector  with  one  element  for  each  vehicle. 
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Its  use  is  as  a pointer  to  ICAR  (i.e.,  vehicle  J 
is  initially  in  route  segment  ICAR(SR0W1 (J) , J) . ) • 

4.  SR0W2  - a vector  with  one  element  for  each  vehicle. 
Its  use  is  as  a pointer  to  PFILE  (i.e.,  vehicle  J 

is  intially  at  mission  profile  point  PFI LE (SR0W2 (J) , 
1 to  6)  ) . 

5.  TYM  - a vector  with  one  element  for  each  vehicle. 

Its  use  is  to  record  remaining  station  dwell  times 
for  the  vehicles.  It  is  initialized  to  0 usually. 

This  completes  the  input  data  preparation.  The  program 
operates  with  the  user  providing  interaction  at  a terminal 
(FORTPJVN  device  TTY) . The  user  is  prompted  to  type  in  at  the 
terminal  the  following  file  names: 


ROW  1 
ROW  20 
ROW  34 

ROW  407 


ROUTE  SEGMENT  1 


ROUTE  SEGMENT  20 


ROUTE  SEGMENT  34 


ROUTE  SEGMENT  407 


FIGURE  6-3  A COLLECTION  OF  ROUTE  SEGMENTS 
IN  ARRAY  PFILE 
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b) 


ICAR 

COLUMN  J 

0 

20 

* 

34 

1 

0 

FIGURE  6-4  THE  SEQUENCE  OF  ROUTES  FOR  VEHICLE  NUMBER  j 
GRAPHICALLY  IN  (a)  AND  AS  TABLE  ICAR  IN  (b) . 
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1. 


a file  name  for  a disk  file  containing  the  guide- 
way data. 

2.  a file  name  for  a disk  file  containing  the  mission 
profile  data. 

Subsequent  prompting  occurs  as  the  user  is  asked  to  supply 
a time  interval  before  the  next  flow  calculations.  Enter- 
ing a negative  time  terminates  the  simulation. 

All  output  is  written  on  FORTRAN  device  number  21. 
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7.  CONCLUSIONS 


The  model  described  here  is  capable  of  simulating  a 
fleet  of  electric  vehicles  moving  over  a guideway  network. 
The  output  is  a discrete  and  vehicle/power  rail  interface. 

It  has  been  in  use  for  almost  one  year.  It  is  a useful 
tool  for  predicting  peak  power  demand,  voltage  droop,  exces- 
sive harmonic  currents. 
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Appendix  A 


A Program  Listing  (AUGUST  1979) 


FILENAME 

CONTENTS 

MAINZ. F4 

the  main  program 

PARANR.F4 

the  common  block  and  declarations  used  in  each 

module 


SETUP. F4 

subroutines  SETPEL  AND  SETGW 

VEHDYN.F4 

subroutine  VEHDIN 

MOTOR. F4 

subroutine  MOTOR 

TEFF.F4 

block  data 

NR2.F4 

subroutine  NR 

F0RMY.F4 

subroutine  FORMV 

NEWPQ.F4 

subroutine  NEWPQ 

PCU.F4 

subroutine  PCU 

GPJ.F4 

subroutine  GETJ,PUTJ 

DISTOR.F4 

subroutine  DISTOR 

OUTPUT. F 4 

subroutines  RITE  § VEHSUM 

HEDWAY.F4 

subroutine  HEDWAY 

NUPNT.F4 

subroutine  NUPNT 

NUGWQ 

subroutine  NUGWQ 
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Appendix  B 

A Program  Concordance 

The  following  table  lists  the  important  terms  and  variables 
used  in  the  program  which  is  listed  in  Appendix  A: 

1.  Main  program  - it  prompts  the  user  for  input  data,  sets  up  the 
state  of  the  vehicles,  initiates  the  load  flow  computation  and 
updates  the  state  of  the  vehicles.  The  latter  two  steps  can 
be  repeated.  Main  global  symbols: 

WEL(II)  - vehicle  II  velocity  (MPH) 

P0S(1,II)  - location  of  vehicle  II  in  terms  of  a guideway  segment 
identifier 

P0S(2,II)  - vehicle  II  position  (FEET)  from  a reference 
TYM(II)  - vehicle  II  dwell  time  (sec.) 

VACC(II)  - vehicle  II  acceleration  (MPH/sec) 

TVM(II)  - motor  terminal  voltage  at  vehicle  II  (VOLTS) 

TIM(II)  - motor  terminal  current  at  vehicle  II  (AMP) 

TTM  - if  POWER  FLOW  INTO  MOTOR  THEN  0,  ELSE  NONZERO. 

2.  SETPFL  - a subroutine  which  prompts  the  user  for  a source 

(disk  file  name)  of  the  mission  profile  input  data. 

IFIL  - mission  profile  data  file  name 

3.  SETGW  - a subroutine  which  prompts  the  user  for  a source  (disk 

file  name)  of  the  guideway  input  data. 

IFIL  - guideway  data  file  name 

4.  VEHDYN  - a subroutine  which  handles  the  vehicle  dynamics 
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5.  MOTOR  - a subroutine  which  models  a separately  excited  dc  motor 

6.  TEFF  - the  block  data  to  initialize  values  in  common 

ICARS(II)  - number  of  cars  coupled  together  for  vehicle  II 
NUMPAS(II)  number  of  passengers  on  vehicle  II 
JERKM  - Maximum  jerk  limit  (MPH/s^) 

ACCM  - maximum  acceleration  limit  (MPH/s) 

7.  NR  - the  Newton-Raphson  load  flow  subroutine.  It  will  iterate 

5 times  or  until  the  power  mismatch  is  less  than  EPSIL 

8.  FORMY  - a subroutine  called  by  NR  to  calculate  the  admittance 

matrix  Y 

9.  NEWPQ  - a subroutine  called  by  NR  to  calculate  the  scheduled  P 

and  Q values  for  each  vehicle. 

GLOBALS 

P(II)  - fundamental  real  power  (per  unit)  at  vehicle  II. 
Q(II)  - fundamental  reactive  power  (per  unit)  at  vehicle 
II 

10.  PCU  - a subroutine  called  by  NEWPQ  to  model  the  power  condition- 
ing unit  which  converts  3 phase  power  into  dc  power 
for  the  motor 
LOCAL 

PRUX  - constant  auxiliary  load  real  power  at  each 
vehicle  (watts) 

QUAX  - constant  auxiliary  load  reactive  power  at  each 
vehicle  (volt-amp) 

AMU  - constant  commutation  angle  (radians) 
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GLOBAL 


ALPHA  - SCR  firing  angle  (radians) 

BP  - fundamental  real  power  (per  unit)  at  vehicle  under 
consideration 

BQ  - fundamental  reactive  power  (per  unit)  at  vehicle 
under  consideration 

11.  GETJ (VAL , I I , J , K)  - a subroutine  used  by  NR  to  look  up  Jacobian 

values.  It  assigns  to  VAL  the  Jacobian 
element  at  row  II,  column  J.  If  K is  zero 
then  the  search  for  the  Jacobian  element 
begins  at  the  end  of  row.  Otherwise  the 
search  begins  at  JCOB(K). 

12.  PUTJ (VAL , I I , J ,K)  - a subroutine  used  by  NR  to  update  Jacobian 

values.  It  assigns  to  the  Jacobian  element 
at  row  II,  column  J the  value  VAL.  If  K 
is  nonzero,  then  begin  the  search  for 
Jacobian  element  at  row  II,  column  J at 
JCOB(JNC(K))  - otherwise  at  JCOB ( JRH (K) ) . 

13.  DIST0R  - a subroutine  which  evaluates  the  Fourier  coefficients 

of  the  harmonic  currents.  Distortion  currents  at  the 
vehicle/power  rail  interface  and  at  the  electric 
utility/network  interface  are  calculated. 

LOCA-NHAR(I)  - a table  of  harmonics  to  consider 
GLOBAL 

DG(I)  - RMS  harmonic  current  (per  unit)  at  node  I. 

14.  RITE  - a subroutine  which  prints  out  the  power  flow  summary 

15.  VEHSUM  - a subroutine  which  prints  out  the  state  (of  the 

vehicles)  summary 
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16. 


HEDWAY  - a subroutine  which  calculates  the  distance  between 


vehicles.  Its  output  is  a message  that  the  mission 
profile  and  vehicle  performance  have  caused  some 
intervehicle  distances  to  shrink  below  a limit 
(called  HDWY) . 

LOCAL 

HDWY  - the  limiting  distance  (feet)  in  front  of  a 
vehicle.  It  is  the  vehicle  length  (feet) 
plus  vehicle  velocity  (feet  per  sec)  times  a 
reaction  time  (sec)  plus  emergency  decelera- 
tion ( feet/sec/sec)  times  velocity  squared 
divided  by  2. 

TAU  - vehicle  reaction  time  (sec) 

EACCEL  - emergency  deceleration  rate  (feet/sec/sec) 

VLNGTH  - vehicle  length  (feet) 

17.  NUPNT  - a subroutine  which  advances  the  vehicle  fleet  to  a 

new  point  on  the  mission  profile 
GLOBAL 

DT  - amount  of  time  to  advance  (sec) 

18.  NUGWQ  - a subroutine  which  updates  all  the  vehicles’  positions 

on  the  guideway  segments.  It  changes  the  vehicle 
numbers  stored  in  array  GWAY  columns  6 to  10. 
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The  Fourier  analysis  integrals  described  in  Section  3.6 


yield  the  following  expressions: 


f 
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A(N)- 


I 


COS  g 

cos  a - cos ( a+uT 


- cos  n(a+y-y/2)  ^ cos  n(a-y/2) 

n n 


cos(y/2) 

cos  a - c 0 s ( a+  y ) 


- cos(n-l (a+y-Y/2)  cos ( n+1 ) ( a+y- y/ 2 ) 
2(n-l  ) 2(n  + l ) 


+ cps(T/2,) 

cos  a - cos(a+y) 


cos { n- 1 ) ( a- y/2 ) - cos ( n+1 ) ( a- y/ 2 ) 

2(n-l  ) 2(n  + l ) 


+ 


s i n ( y/  2 ) 

cos  a - CO s ( a+y ) 


sin(l-n)(a+y-Y/2)  sin(l+n)(a+y-Y/2) 

2(l-n)  2(l+n) 


+ sin(Y/2) 

cos  - cos(  + ) 


s i n ( 1 - n ) ( g-  Y/  2 ) 
2(l-n) 


- cos  n(g+y+Y/2)  ^ COS (n ( g+y- y/2 ) 

n n 


- cos  g 

cos  g - cos ( g+y) 


- cos  n(g+y+Y/2)  ^ cos  n(a+Y/2) 

n n 


+ cos(-y/2) 

cos g - COS ( g+y) 


COS ( n- 1 ) ( g+y+Y/ 2 ) 

2TK^i1 


cos|n+1 j (g+y+Y/2) 


- cos(“y/2) 

cos  g - cos ( g+y ) 


cos(n-1 ) (g+Y/2) 
l(n-l) 


cos ( n+1 


I (g+Y/2) 


2(n  + l 


+ _s_in(-Y/2) 

cos  g - cos{g+y) 


si n(l-n) (g+y+Y/2)  s i n ( 1 +n ) ( g+y+Y/ 2 ) 
2{l-n)  2(l+n) 


+ sin (-y/2) 

cos  g - C0S(g+y) 


sin(l-n)(g+Y/2)  s i n ( 1 +n ) ( g+Y/ 2 ) 

2(l+n)  
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Id 

b(n: 

1 = 

7T 

cos 

a 

cos 

a - 

COS 

( a+y) 

cos 

T/2 

cos 

a - 

cos 

( ot+y) 

cos 

y/2 

cos 

a - 

cos 

(a+y) 

s i n 

Y/2 

cos 

a - 

cos 

(a+y) 

sin 

T/2 

cos 

a - 

cos 

(a+y) 

si  n 

n(a+u+Y/2) 

n 

cos 

a 

cos 

a - 

cos 

(a+y) 

cos 

Y/2 

cos 

a - 

cos 

(a+y) 

cos 

y/2 

cos 

a - 

cos 

( ot+y) 

cos 

y/2 

cos 

a - 

cos 

(ot+y) 

s i n 

yZ2 

sin  n( a+y- y/2) 

sin  n ( a- y/ 2 ) 

n 

n 

- 

s i n ( 1 -n ) ( a+y- y/2 ) 

, s i n ( 1 +n ) ( a+y- y/ 2 ) 

2(l-n) 

2(l+n) 

s i n ( 1 - n ) ( a- y/2 ) 

sin(l+n)(a-Y/2) 

2(l-n) 

2(l+n) 

- 

CO  s ( 1 - n ) ( a+y- y/ 2 ) 

cos( 1+n) ( a+y-Y/2) 

2(l-n) 

2(1+n) 

COS ( 1 -n ) ( a- y/ 2 ) 

cos ( 1 +n ) ( a- 

y/2) 

— 

2(l-n) 

2(l+n) 

_ 

sin  n( a+y- y/2) 
n 

sin  n(a+y+Y/2) 

sin  n( a+Y/2) 

n 

n 

s i n ( 1 - n ) ( a+y+Y/ 2 ) 

. s i n ( 1 +n ) ( a+y+Y/ 2 ) 

2(l-n) 

2(l+n) 

. sin(l-n) (a+Y/2) 

. si n ( 1 +n ) (a+Y/2 ) 

2(l-n) 

2(l+n) 

COS ( 1 -n ) ( a+y+Y/ 2 ) 

cos ( 1 +n ) ( a+y+ y/ 2 ) 

2(l-n) 

2(l+n) 

cos ( 1 - n ) ( a+Y/ 2 ) cos ( 1 +n ) ( a+y/2 ) 

2(1 -n)  2(l+n) 
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